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Abstract 

We have calculated the decay r — > vk{K)^. We present the photon energy spectrum, the 
meson-photon invariant mass spectrum and the integrated rate as a function of a photon 
energy cut or an invariant mass cut. Both the internal bremsstrahlung and the structure 
dependent radiation have been taken into account. To this aim we have parametrized the 
form factors Fy and Fa, which determine the structure dependent radiation. Observables 
especially suited for the measurement of the structure dependent form factors are found 
and implications on the width of the a\ discussed. 



1 Introduction 



Present and future high luminosity electron positron colliders and especially tau-charm factories will 
produce very large samples of tau pairs in the near future. This will give the opportunity to observe 
rare tau decays, for example radiative semihadronic modes, the simplest of which are r — ► wiry and 
r — ► uKj, which we will therefore calculate in this paper. 

The decay r — ► wiry is related by crossing symmetry to the radiative pion decay (it — ► ewj), 
where the electron plays the role of the tau lepton. There is, however, a crucial difference, namely 
that the momentum transfer squared t between the leptonic and the pion-photon system, which is 
almost zero in the case of the pion decay, can take any value up to rr? T in the radiative tau decay. 
Several resonances, which can contribute only as small virtual corrections in the pion decay, can 
now be produced as real particles and duly have to be taken into account. 

The decay amplitude can be divided into internal bremsstrahlung, where the pion behaves like 
a pointlike particle and which therefore is completely determined by quantum electrodynamics, and 
the structure dependent part, which is governed by strong interaction effects and is parametrized by 
two formfactors Fy(t) and i*A(t). So the main task is to obtain these formfactors. We will assume 
them to be dominated by several Breit Wigner type resonances and determine their parameters by 
employing different constraints, both theoretical and experimental ones. For instance normalizations 
at t = are fixed by the Adler-Bell anomaly (Fy) and the tt — > ewy measurement (Fa), respectively. 
On the other hand the relative strengths of the different resonances are taken from a high energy 
QCD prediction and experimental data on the formfactors at finite values of t. 

This paper is organized as follows: In Sec. 2 we shortly discuss the kinematics, the differential 
decay rate and lay down our conventions. In Sec. 3 the general form of the matrix element, photon 
spectrum and the decay rate of t — > wiry are calculated. Their dependence on the two formfactors 
Fy(t) and F^it) are worked out. In Sec. 4 these formfactors are parametrized along the lines indicated 
above. In Sec. 5 the numerical results for the photon and invariant mass spectra and the decay rate 
are presented. In Sec. 6 we estimate the decay r — > vKy and in Sec. 7 we discuss our results. 



2 Kinematics and the Decay Rate 

For definiteness and without loss of generality we consider the decay of the negatively charged tau 
into a pion, a neutrino and a photon 

t~(s) — ► v T {q)n~ (p)j[k) (1) 

only and transcribe the formulae to the case of the decay into a kaon in Sec. 7. The kinematics of 
this decay is equivalent to that of the radiative pion decay We can write the differential decay 
rate for an unpolarized tau in such a way that it depends on two Lorentz scalars, taken to be x and 

y- 

2s -k 2s -p 
x := — y := r (2) 

In the tau rest frame x and y are the energies E-y and E n of the photon and the pion, respectively, 
expressed in units of m T /2: 

E i = ~y x E * = ~yv (3) 

The kinematical boundaries for x and y are given by 

< x < 1-r 2 



2 



^.2 

1 - x + < y < 1 + r 2 (4) 

1 — x 



where 
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using the central value of the new r mass experiments [§] 

m T = 1777 MeV. (6) 
Let us quote some useful formulae which relate different kinematical expressions: 
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p-k = ^-{x + y-l-r 2 ) (7) 

t ■= (s - qf = (k+p) 2 = m\(x + y - 1) (8) 

Now we are in position to present the formula for the decay rate of an unpolarized tau decaying into 
|^ r 7T7) The differential decay rate is given by 

where the bar over the matrix element denotes summing over the photon polarization and neutrino 
spin and averaging over the tau spin. 

Choice of the tau rest frame, integration over the neutrino momentum and the remaining angles 
and introduction of the dimensionless variables x and y lead to 

d?T m 



T -AM\ 2 (10) 



dx dy 2567T 3 
The integration over y yields the photon spectrum 

dx Ji- x+ -lL. dxdy 

Because of the infrared divergence of the internal bremsstrahlung a low energy cut must be intro- 
duced for the photon energy, eg. by requiring 

x > x (12) 

The integrated decay rate is then obtained by 

rl—r 2 fir 

r(x ) = T(E ) = J x dx— (13) 

and does depend on the cut in the photon energy (E = —^-xq). 

Instead of the invariants x and y we will also use x and z where z is the scaled momentum 
transfer squared: 

z = -^r = x + y-1 (14) 
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Their kinematical boundaries are given by 

„2 



z — r 2 < x < 1 



r 



r z < z < 1 (15) 

Integration of d 2 T/dx dy over x then yields the spectrum in z, ie. the spectrum in the invariant mass 
of the pion-photon system: 

~r(z) = -ri^i) = / dx——(x,y = z-x + l) 16) 

az az Jz-r 2 dx dy 

The integrated rate for events with t > to is then given by 

r 1 dF 

r(z ) = T(VT ) = J dz—(z) (17) 



where, of course, zq = — Note that the cut in t acts as both an infrared and a collinear cut off. 

3 General Structure of the Matrix Element and the Decay Rate 

In complete analogy to the case of the radiative pion decay Q, the matrix element for the decay 
r — * vttj can be written as the sum of four contributions: 



with 



m[t (s) —>Mq)-x~(ph(k)] =m IBt +m IBw +Mv + Ma (is) 



Mib t = -G F cos6 c ef 1T p iJ ,e v (k)L^ 
Mib„ = G F co S 9 c eUe v (k) f 2 P^+P^ _ g ^ ) jy 

My = ^^^ie^(k)k^ 



G F cos 6 c eF A (t) 
\/2 m„ 



M A = * ° ° r J e v (k) ((p„fc„ - (p ■ k)g„ v ) L» (19) 



where G F is the Fermi constant, 9c the Cabibbo angle, e the electric charge and e v the polarisation 
vector of the photon. For our numerical results we use the pion coupling constant f w to be 0, || 

/ w = 92.5 MeV (20) 

Fy and Fa are the so called structure dependent form factors which are parametrized in the next 
section. Note that our definition implies that these form factors are dimensionless. Finally and 
If' 1 ' are leptonic currents defined by 

"~ = ^>^- (L~^-V " (s) <21> 

In fact the four terms correspond to the Feynman diagrams in Fig. 1: 
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• Mj Bt i s the bremsstrahlung off the tau (a), 

• MiB n the sum of the pion bremsstrahlung (b) and the seagull (c), 

• My is the structure dependent vector (e) and 

• Ma the structure dependent axial contribution (f). 

Note that all ambiguities due to the strong interactions have been parametrized in terms of the two 
form factors Fa and Fy. In fact these form factors are the same functions of the momentum transfer 
t as those in the radiative pion decay, the only difference being that t now varies from up to m T 
rather than just up to m n . 

The two matrix elements Mi Bt an d Mi BlT are not separately gauge invariant, but their sum, 
ie. the (total) matrix element for internal bremsstrahlung IB 

M ib = M IBt +M IBw (22) 

is, as are M.y and Ma- We also define the (total) structure dependent radiation SD, viz. by 

Msd = M v + Ma (23) 

Using standard technics the matrix elements may be simplified, the main rearrangement being 

We obtain 



Mib = G F cos 9 c ef 1 ,m T u{q)^ + 



P- e + W s ■ e 



u(s) 



p ■ k 2s ■ k s ■ k 

Msd = ° F °% L^L^k^^- + u{q) 1+ [(p • k)4 - (e • p)k/] u(s)^\ (25) 
The square of the matrix element is then given by 



\M\ 2 = \M IB \ 2 + 2Re(M IB M* SD ) + \M S d\ 2 , (26) 

where — as stated in the previous section — the bar denotes summing over the photon polarization 
and neutrino spin and averaging over the tau spin. 

In order to make our results transparent we have divided the decay rate into several parts: 



the internal bremsstrahlung part Ti B arising from \Mi B \ 2 , 



• the structure dependent part Tsd arising from |A^5d| 2 



• the interference part Tint arising from 2Re(Mi B M* SD ). 

Furthermore we have subdivided in the obvious way into Tyy arising from the vector matrix 
element squared, Taa form the axial part and the vector-axial interference term TyA and similarly 
T/att gets divided into the internal bremsstrahlung- vector interference Ti B _y and the internal 
bremsstrahlung-axial interference Ti B ^a- 
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So our master formulae are: 



^total 

One finds for the differential decay rate 



+ Tsd + ^INT 
Tyy + TyA + Faa 
Fib-V + Fib-A 



(27) 



d 2 T IB 


a 






dx dy 


,r 


>{1 -r 2 ) 2 


d 2 T SD 


4 


'\F v \ 2 f vv (x, 


dx dy 


16vr /2 m 2 


d 2 FiNT 


a m 2 




fiB-v(x,y 


dx dy 




where 









\F v \ 2 f vv (x, y, r 2 ) + 2Re(F v FX)fvA(x, y, r 2 ) + \F A \ 2 f AA (x, y, r 2 



fiB-v(x, y, r 2 )Re(F v ) + f IB -A(x, y, r 2 )Re(F A ) 



(1-r 2 ) 2 



(1 -r 2 ) 2 
(28) 



, , 2x [r A {x + 2) -2r 2 (x + y) + (x + y -l){2 -2,x + x 2 + xy)]{r 2 -y+l) 

JiB{x,y,r ) = , i\2 2 

(r z — x — y + l) z x^ 

fw(x,y,r 2 ) = -[r\x + y) + 2r 2 (l-y)(x + y) + (x + y-l)(-x + x 2 - y + y 2 )] 
f AA (x,y,r 2 ) = f vv {x,y,r 2 ) 

f VA (x,y,r 2 ) = [r 2 (x + y) + {I - x - y){y - x)){r 2 -x-y + l) 



fiB-v(x,y,r 2 ) = - 
f IB -A(x,y,r 2 ) = - 



(r 2 - x - y + \){r 2 - y + I) 

x 

4 - 2r 2 (x + y) + {l-x + y){x + y-l)]{r 2 -y + l) 
(r 2 — x — y + l)x 



(29) 



In the approximation r 2 rs (vanishing pion mass) the formulae simplify to 



f IB (x,y,0) 

fvv(x,y,0) 
f VA (x,y,0) 

fiB-v{x,y,0) 
fiB-A{x,y,0) 



[l + (l-x) 2 -x(l-j/)](l-j/) 

[x + y — l)x 2 
(x - x 2 + y - y 2 ){x + y - 1) 
(x + y- l) 2 (y-x) 
(x + y-l)(l-y) 

x 

(x-y-l)(l-y) 

x 



(30) 



Note that we expressed the radiative decay rate in terms of the rate of the non-radiative decay 
(r -> ^ T vr)): 

„2\2 



Gf.cos 2 ^/ 2 
8vr 



m%{l 



(31) 



Although the factor m 2 /(/ 7r m 7r ) in front of Fa and Fy is rather large a number (« 240), we should 
mind that the form factors F A and Fy themselves are of the order of 10~ 2 (see next section), so 
that we expect the different contributions to be of the same order of magnitude. Note also that 
the internal bremsstrahlung and the interference term are infrared divergent in the limit x — > 
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(vanishing photon energy). Furthermore there is an enhancement of the internal bremsstrahlung 
(but not of the interference term) in the region 



(x + y - 1) — » r 2 44> t 



(32) 



This enhancement is due to large logs of the pion mass, the usual collinear divergencies for massless 
particles. 

Next we would like to make a short comment on the internal bremsstrahlung part of the decay 
rate. Let us first give the corresponding expression in the radiative pion decay || (neglecting m e ) 



d 2 TiB{Tt — ► ewy) a 
dx dy 2ir 



r(7r — ► eu)IB(x,y) 



with 



IB(x,y) 



(l-y)[l + (l-x) 2 ] 



x 2 {x + y — 1) 

where now, however, x and y are the photon and electron energies in units of m n /2: 



-x 



(33) 
(34) 

(35) 



Expressed in this form eq.(34) is identical with eq. (|30|) , apart from the x(l — y) term. On the 
other hand keeping in mind the definition of x we note that m T /2 and m T /2, respectively, set the 
scale for photons to be considered as "hard" or "soft". This means that the formulae for internal 
bremsstrahlung should be similar for radiative tau and pion decay, once they are expressed in terms 
of x and y, which is indeed the case. So photons from the r- or n- decays of comparable "softness", 
e.g. comparable x, have very different energies. 

Finally we present analytical expressions for the invariant mass spectrum: 



dT IB 

dz 



a 
2^ 



[l-z) + 2r 2 (z - z 2 ) -Az + 5z 2 - z 3 + 

i r T 



+(r 4 z + 2r z z - 2z - 2z z + z 6 ) In z 



dTyy 


a m\ (z — 


\) 2 {z-r 2 


) 3 (l + 2z) 


dz 


48vr /2 m 2 


z 2 




dT V A 


= 






dz 






dT A A 


a m\ (z — 


\) 2 {z-r 2 


) 3 (l + 2z) 


dz 


48vr /2 m 2 


z 2 




dTiB-v 


a m 2 (z 


-r 2 ) 2 (l- 


z + z\nz) 



z 2 — r 2 z (1 — r 2 ) 2 
\Fy\ 2 r ™ 



(1 - r 2 ) 2 



dz 
dTiB-A 



F A \ 



(i 



„2\2 



2y/2lT fnm^ 



a mz 



Re(Fy) 



(l-r 2 )2 



dz 



2y/2n Urn* 
+ (r 2 z — 2z — z 2 ) In z 



r 2 (l - z) - 1 - z + 2z 2 + 



z — r 



(1 



-2\2 



(36) 



Note that the interference terms IB-V and IB- A are now finite in the limit z — > r 2 , which proves 
that their infrared divergencies are integrable. 
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4 Parametrization of the Form Factors 



4.1 The Axial Form Factor 

The axial form factor F A (t) fulfills a once subtracted dispersion relation. Since this form factor is 
measured at t = in the radiative pion decay we will use this value to fix the subtraction constant. 
Furthermore it is supposed to be dominated by the a\ meson. The following form includes these 
constraints: 

F A (t) = F A (0)Res ai (t) (37) 
where F A (0) is taken from The Particle Data Group § 

F A (0) exp = -0.0116 ± 0.0016. (38) 

and Res ai (t) is a normalized Breit Wigner resonance factor 

Res ai (t = 0) = 1 (39) 

Note that most experiments on the radiative pion decay have been performed with pions at rest. 
These experiments are sensitive mainly to (1 + 7) 2 only, where 

7 = ^ (40) 

So they also allow for a positive sign solution for Fa(0: 

F A (0) = +0.0422 (41) 



Measurements performed with pions in flight support our value in (pq). However, we will also 
comment on implications of the other solution. 

Since the decay of the a\ is strongly dominated by its three pion decay we approximate in the 
Breit Wigner the total width by the three pion width. Furthermore we will use a energy dependent 
width as given in Q 

m 2 

mi -t- im ai T ks (t) 



where 



with 



r ks (t) = -j%r ai (43) 



if t < 9ml 

g(t) = { 4.1(t-9m 2 r ) 3 [l-3.3(t-9m 2 r ) + 5.8(t-9m 2 r ) 2 ] if 9m 2 < t < {m p + m w ) 2 \ (44) 
t(1.623 + 10.38/t 2 - 9.32/t 4 + 0.65/t 6 ) else 

(all numbers in appropriate powers of GeV). 

We have studied the influence of a energy-dependent width by comparing our results to those 
obtained from 



m 2 - im ai T ai 

n 

with a constant width T, 



Res ai (t) = 2 ai . 1 I (45) 
m l -t- im ai T ai 



ai • 
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In order to check our form factor we use it to compute the radiative decay of the a\. Indeed in 
the narrow width approximation, ^(m^ ) can be related to r ai ^ 7r7 , once the transition amplitude 
for r — * va\ is known. Let this transition be parametrized by the coupling constant f ai 

M{r vax) = G F cos 6 c f ai e^ (46) 

Then 

Using equation (j37]) , this can be reexpressed in terms of Fa (0) as 

r — 7 = ^-3^l^(0)| 2 fl-^V (48) 



48 mlfl 



m 



From the experimental result with || 

{r ai ^j) exp = 640 ± 246 keV (49) 

we obtain for the central values 

f ai = 0.09 GeV 2 (50) 

which has to be compared with 

f ai =f p = V2m p f n = 0.10 GeV 2 (51) 

from chiral symmetry and the Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin (KSRF) relation (8| 
or with f ai extracted from r — > a\v data: 

,2 m2 ai m2 r 1 BR{t -> OiU) 

/ai 24vr 2 cos 2 %^_^2x2^ + 2 ^^ 2 x BR(t ^ e~uu) 1 ) 



"It 



(where BR denotes branching ratio). With present data, BR(t — ► va\) = 6.8% || we obtain 

f ai = 0.131 GeV 2 (53) 

from the direct measurement. 

So we observe a reasonable agreement between the predicted and measured radiative decay width 
of the a\. 

4.2 The Vector Form Factor 

The hypothesis of the conserved vector current (CVC) relates the vector formfactor Fy to the 
formfactor i ? 7r7 of the vertex 7r°77* by 

Fy{t) = - m - F J2 {t) (54) 

where F ny (t) is defined by 

M [ir°(p) — 7 (fc)7*(?)] = ie^^e^V^)^) (55) 

The case q 2 = corresponds to the decay of a pion into two real photons. 
There are three constraints for the form factor F n ^: 
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there is a low energy theorem for i ? 7I - 7 (0), 
the slope of F n ^(t) at t = has been measured and 
there is a perturbative QCD theorem for q 2 - 



oo. 



In good agreement with the measured rate for the decay 7r° 
Adler-Bell-Jackiw anomaly pOl as 



1 



77, -^7(0) is predicted by the 



(56) 



This yields 



F v (0) = -0.0270, 



(57) 



which is almost within one standard deviation of the experimental value || measured in the radiative 
pion decay (where t ~ 0), which in our phase convention reads 



F V (0) 



Perturbative QCD predicts 



exp 



F n Jq 2 



-0.017 ±0.008 



-00 



2/tt 
^2 



(58) 



(59) 



Furthermore the slope of F^~ at t = has been measured. Using the definition 



F^{t) = ^(0) 



1 + ci- 



rri 



+ 0(t 2 



the experimental value || is 



a exp = +0.0326 ± 0.0026. 



(60) 



(61) 



Assuming to be dominated by p like resonances, we will try to parameterize F^{t) in 
accordance with the constraints described above using the following models: 



monopol 



dipol 



tripol 



Fir^y (t) 



Resx {t) 



+ A 



? 7r 7 (i) = 47i _ 2 1 [Resp(t) + ARes p /(i) + //ReSp»(t)] y 



1 



+ A + ,u 



(62) 
(63) 
(64) 



All these models automatically satisfy the low energy constraint, because the resonance factors are 
again normalized according to 

Resx(0) = 1 (65) 



Furthermore we suppose 



Resx (t — ► — 00) 



m 



x 



(66) 
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This is fulfilled by 

Roo„W — 



2 

Res x (t) = —2 m , x - , . (67) 



m x ~ t ~ irnxT pw (t) 

with an energy dependent width corresponding to the P-wave two body phase space 

if t < Ami 

3 



Fpwit) = { m x \Jt ~ Am\ 



The form 

ml-im x T x 




m\ — Ami. 



Tx else 



m x — t — imx^x 

with a fixed width Tx will again be considered for comparison. 

If we use the physical rho mass m p = 768 MeV in the monopol parametrization, the high energy 
limit is 14 % below the QCD prediction. Brodsky and Lepage note in 12] that the choice 



m x 



M BL ■= 2v / 2?r/ 7r = 822 MeV (70) 



satisfies the low and the high energy constraint simultaneously with a monopole formfactor. The 
slope parameter a depends on the resonance mass as follows: 

fl - m 2 / 2 _/ 0.0309 if m x = m p 

«-"V/m x -j 02?0 if mx = M BL (71) 
So mx = Mbl is not a very good approximation in the low energy region and therefore we will use 

ni x = m p (72) 

for the resonance factor of the monopole parametrization. This monopole version defines the model 
N =1. 

In the dipol parametrization the parameter A must be 

A = = 0.0584 (73) 

m 2 p ,-M BL 

in order to satisfy the QCD theorem. (If not stated otherwise, particle parameters of || are used 
everywhere.) For the slope parameter this implies 

a = TXT ( —2 + A) = °- 0296 ( 74 ) 
1 + A \m 2 p mp J 

in satisfactory agreement with (|6l|). The dipol formfactor defines the model N = 2. 

In order to determine the parameters A and (i of the tripol model, two more constraints will be 
considered, viz. the experimental values for r p ^ 7r7 and T p >^ n0J . Using the narrow width approxima- 
tion, parametrizing 

M (t - pv) = G F cos ej^L" (75) 



and dividing Fy into 



F v (t) = + F^ p,) + Ff ] (76) 
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in an obvious notation, T p ^^ can be calculated in terms of the parameters of Fy\ The result is 



r 



a m' 



inz 



1536vr4/2 /p 2 \ ml (1 + A + p) 



2 ' 



We will use f p as given by the KSRF relation, 

f p = V2m p f w = 0.10 GeV 2 
Equally well we could extract f p directly from t ^ pv data: 



r 



1 



24vr 2 cos 2 6, 



C 



( m p\ 



BR(t -> pv) 



With present data f| this is 



/„ = 0.11 GeV 



We have checked that the use of this value would not change the final result significantly. 
So for A = p = (monopole) we obtain 

(r p^irj) monopole — 80.5 keV 

On the other hand the experimental number is 

(r p ^) exp = (68 ± 7) keV 

and so 

l + X + p-- 



(-Tp— >7T7 ) monopole 



(r 



p— fiT'y Jexp 



1.088 



,(p') 



(77) 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



In order to relate IV —►t^ to the parameters of Fy , we assume Fy to be given by the Feynman 
diagram of Fig. 2. So we need f p >, g p 'ujn and / w . / p > can be extracted from a fit to e + e~ — > 27r data 
0], which may proceed through p, p' or p" resonance channels (see Fig. 3). Neglecting the pion mass, 
g pn7T and g p '^ are related by 

dp'lTTT 



gp-n 



/ Fp'^27r m p 

r p ^27r ny 



0.23 



where we used Ipfl 



Fp'^ 2 7r = 0.24 x 0.21 x 310 MeV 



Let the fit to e + e — ► 2ir be parametrized by 



(ReSp(t) + crRes„/(t) + pRes p n(t)) 

l + O" 



+ P 



Then we have 



m 2 : (crm p ,) = (f p g p7r n) ■ (/p'V^) 
Reference gives (see their model N = 2, table 1) 



(84) 
(85) 

(86) 
(87) 



a = -0.103 
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Therefore 







9p7V7T 


fp 




9p f 7T7T 



ml 



a = 1.63 



The coupling g p / uw is defined by 

M(p'{q) -> u{k)ir{p)) = v w e F /(?)£ , '(i)^ 
and can be deduced from 



bp 1 



f^ m 3 i _ ^ 

V m P '. 



96vr 



0.21 x 310 MeV 



as 



Assuming = / p /3, comparison of the diagram of Fig. 2 and Fy then gives 



|A| 



2^2 



1 + X + fj.) 



m>2 Fy(0) 



/p 



0.136 



(90) 

(91) 
(92) 

(93) 



If A and 1 + A + /i are given with moderate errors, \x can only be calculated with a large uncertainty. 
So instead we use the QCD limit to express // in terms of A and 1 + A + [i as 



1 



I* 



m. 



[(1 + A + ix)Ml L - ml - \m\ 



We find 



A = +0.136 [i = -0.051 =► 1 + [i + A = 1.085 
A = -0.136^ /x = 0.15 => l + /i + A = 1.014 



(94) 



(95) 



So in order be consistent with 1+//+A = 1.088 from F„-+ n j, A must be +0.136. The slope parameter 
a then is 



mt 



1 



A 



1 + A + \i \m p + m 2 pl + m 2 n „ 



0.0291 



(96) 



in satisfactory agreement with equation (|6l] 
So for the tripol model, 



A 
/' 



+0.136 
-0.051 



(97) 



is a choice for the two parameters which is consistent with all the four constraints (the QCD limit, 



the slope a, T p -^~ and F 



p — fKijJ , 



and defines the model N = 3. 



5 Results for the Spectra and the Decay Rate 

We consider as standard the following choices: 

• energy dependent widths in both form factors, 

• model N = 3 (tripol) for Fy, 
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• F A (0) = -0.0116, and 

• T ai = 0.4 GeV. 

In Fig. 4 (a)-(d) the resulting photon spectrum is displayed for 0.1 < x < 1. For "soft photons" 
(xq < 0.2) the internal bremsstrahlung completely dominates (Fig. 4(a)). Note that for very soft 
photons the multi-photon production rate becomes important, and so our order a results are not 
reliable too close to the infrared divergence x = 0. We have checked our internal bremsstrahlung 
spectrum against the data produced from PHOTOS and we obtain excellent agreement, apart from 
the region of very large x. 

The spectrum is significantly enhanced by structure dependent effects for hard photons (xq > 
0.4) (Fig. 4(b)). In Fig. 4(c) we show that the vector resonances dominate the structure dependent 
part while in Fig. 4(d) we plot the interference between bremsstrahlung and structure dependent 
part. 

While the integration over the internal bremsstrahlung needs an infrared cut-off, the structure 
dependent part does not. In order to get a feeling on its importance we have performed the integration 
over the complete phase space which yields 

T vv = 0.75(0.86) • 10~ 3 r r ^ 

v AA = o.33 • io- 3 r T ^ w 

T V A = 

^T SD = 1.08(1.20) • 10- 3 r T _» w (98) 

using the tripol parametrization for Fy. The numbers in parentheses are the predictions of the 
monopole fit. 

With the value of 0] for the branching ratio BR(r — > ttv) this translates into 

BR vv (t i/7T7) = 0.87(1.00) • 10" 4 

BR aa (t -» vwi) = 0.38 • 10 -4 (99) 

These values may be compared to narrow width estimates: taking into account the lowest lying 
resonances p, a\ we obtain 

BR(t -> up -► virj) ps BR(t -> up) x BR(p -► vr 7 ) ps 22% x 4.5 • 10~ 4 = 0.99 • 10 -4 

BR(r -» uai vn-y) w BR(t i/ 0l ) x BR(ai tt 7 ) ps 6.8% x 1.6 • 10~ 3 = 1.01 • 10~ 4 (100) 
Here we have used the central value 

BR(ai - 7T 7 ) = 640 keV = 1.6 • 10- 3 (101) 
v n 400 MeV v ; 

So for both the vector and the axial vector channel the narrow width approximation and the picture 
of a decay chain r — > vp(ai), p(a\) — * vr 7 give reasonable results (do not forget the large error on 
BR(ai -» vr 7 )). 

In tables || and || we display for two values of the photon energy cut how the different parts 
contribute to the total rate (using the standard parameters). Note that the only term involving 
Fy which is of importance is Tyy, while for F A the internal bremsstrahlung-structure dependent 
interference part Tib- a is of the same size as T AA . While the vector-axial interference Ty A can 
always be neglected, the internal bremsstrahlung-structure dependent interference Tint is rather 
small but not negligible. 
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Table 1: Contribution of the different parts to the total rate, using E~ > 400 MeV 



T vv (xq) 
T AA (x ) 

Fib-a(xq) 



0.55 • io~ 6 r T ^ U7: 

0.01 • lO -3 ^^ 

o.3i • io~ 3 r T ^ l/7r 

0.07 • 10~ 3 r T _ ww ' 

0.35 • io- 3 r T _„~ 



^int(xq) 



= 1.48 • kht t 



0.87- 10" 3 r T 



0.42 • io _3 r r ^^ 



2.76 • io _3 r T 



Table 2: Contribution of the different parts to the total rate, using > 50 MeV 



r vv (x ) 

T AA (x ) 

Tib-v(xo) 
Y IB -a{xq) 



= 0.75 • i(r d r T _ 
= o.oo • io _3 r T . 
= 0.33 • io~ 3 r r _ 

= 0.05 • io~ 3 r T ^ 
= o.57-io~ 3 r T ^ 



Tib(xo) 
> Tsd{xo) 

^int(xq) 

Ftotal(xo) 



i3.i • io _3 r T 
i.08 • io~ 3 r T 

0.62 • 10" 3 r T 
14.8 • io _3 r T 
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Table 3: Dependence of the total rate on variation of the parameters. Standard choices are implied 
wherever not stated otherwise. E~ > 400 MeV. 



Variation 


1 total I J-U 1 T^UTT J 


N = 1 (monopol) 


2.88 


N = 2 (dipol) 


2.77 


N = 3 (tripol) 


2.76 (standard) 


fixed widths for the p's 


2.80 


r oi = 0.6 GeV, variable width 


2.60 


r oi = 0.4 GeV, variable width 


2.76 (standard) 


r oi = 0.25 GeV, variable width 


3.05 


r 01 = 0.6 GeV, fixed width 


2.64 


r oi = 0.4 GeV, fixed width 


2.79 


r ai = 0.25 GeV, fixed width 


3.07 


F A (0) = -0.0100 


2.63 


F A (0) = -0.0116 


2.76 (standard) 


F A (0) = -0.0132 


2.90 


F A (0) = +0.0422 


4.94 


F A {t) = F A (0), F v {t) = F v (0) 


1.25 



In table ^ the dependence of the total rate on the variation of the parameters is displayed for 
E~ > 400 MeV. It turns out that the uncertainties in the vector form factor Fy (whether we use a 
monopole, dipole or tripol form and whether we assume fixed or energy dependent width) only have 
rather little impact on the total rate. Its uncertainty is dominated by the ambiguities of the axial 
form factor F A , in which connection the value for the a\ width and the normalization error at t = 
are of comparable importance. 

Note that choice of the positive sign solution for F A almost doubles the integrated rate, and so 
by measuring r — * virj a clear decision for one solution is possible. 

In the last column of this table we display the result of an integration with the formfactors 
fixed at their values at t = without Breit-Wigner resonances. It turns out that the Breit-Wigner 
enhancement substantially increases the total rate. 

Fig. 5 displays the integrated decay rate in variation with the photon energy cut Eq (E^ > Eq). 
We give the prediction using the standard parameter set, a lower and upper limit and the contribution 
from pure internal bremsstrahlung. The lower limit is obtained using 

• N = 3, tripole form of Fy(t), 

• T ai = 0.6 GeV, 

• variable widths for the formfactors, and 

• F A (0) = -0.0100 

the upper limit corresponds to the following choice: 

• N = 1, monopole form of Fy(t), 

• T ai = 0.25 GeV, 
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Table 4: Contributions to the integrated rate T(y/to) using the standard parameter set (all numbers 
for the rate in units of 10~ 3 T T ^ U7T ) 



Vk( MeV) 


Tib 


T V v 






FlB-A 


r total 


400 


3.19 


0.74 


0.33 


0.05 


0.51 


4.83 


700 


0.77 


0.61 


0.33 


0.13 


0.33 


2.17 


1000 


0.19 


0.09 


0.28 


0.07 


0.08 


0.72 



• fixed widths for the formfactors, and 

• F A (0) = -0.0132. 

In Figs. 6 (a) — (d) we show the pion-photon invariant mass spectrum. We find that a much 
better seperation between internal bremsstrahlung and structure dependent effects is obtained here 
(Fig. 6 (a) and (b)), as compared with the photon spectrum. So this spectrum (in the region 
Vi > 700 MeV) is an observable that is much better suited for studying the form factors than 
the photon spectrum is. It turns out that with the a\ width of T ai = 400 MeV, the p is the only 
resonance visible in the invariant mass spectrum. The total structure dependent (SD) spectrum is 
the sum of the vector and axial vector Breit-Wigners (VV and AA), as the vector-axial interference 
(VA) vanishes in the invariant mass spectrum after integration over the other kinematical variable. 
So in the SD spectrum there is a very soft and broad a\ resonance bump (Fig. 6 (c)). The internal 
bremsstrahlung-structure dependent interference radiation (INT) near the a\ is dominated by the 
axial part (IB-A). It changes sign at the a\ mass and rises strongly with decreasing \ft below m ai 
(Fig. 6 (d)). The internal bremsstrahlung (IB) also rises below m ai , and so the a\ bump is completely 
erased in the total spectrum (Fig. 6 (b)). 

In Fig. 7 we plot the spectrum in the a\ mass region using r ai = 250, 400 and 600 MeV. The 
result is that the a\ resonance peak is visible if and only if the a\ width is small (r ai ~ 250 MeV). 

We close this section by discussing the integrated rate with an invariant mass cut. For this 
invariant mass cut vh) we suggest three different values: The rather low value of 400 MeV, which 
gives a rather high integrated rate, the value 700 MeV, which suppresses the internal bremsstrahlung 
without loosing much structure dependent radiation and 1000 MeV in order to focus on the axial 
channel (see Tab. |j). 

For the total rate we find 

r(^ = 400 MeV) = 4.83 (4.55... 5.62) lO^T,-^ 
r(^ = 700 MeV) = 2.17(1.91... 2.95)10~ 3 r T _ +w 
r(v^ = 1000 MeV) = 0.72(0.52... 1.28) lO^T,-^ (102) 

where the central values correspond to the standard parameter set and the lower and upper limit 
are obtained with the same parameter variations as above. 

If we vary T ai from 250 MeV to 600 MeV and keep the other parameters fixed at their standard 
values, the total rate r( v / to = 1000 MeV) varies from 1.00 to 0.56 10~ 3 T T ^ V7r , if we vary F A (0) 
from —0.0100 to —0.0132, the total rate varies from 0.64 to 0.81 10~ 3 T T -> un , So because of the 
experimental uncertainty of F A (0) it is difficult to get much information on the a\ width by measuring 
the integrated rate and the invariant mass spectrum should be used, as discussed above. 
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6 Estimation of the Decay r — > vKj 

Compared with the pionic mode, the decay r — ► zvif 7 is both Cabbibo and phase space suppressed, 
i.e. it is even more rare. So a reasonable estimate for the total rate T^-(xo) will be sufficient. 
For r — * vK'y, in Sees. 2 and 3 the following substitutions must be made: 

cos 9 C — * sin 9 C 

U — /* = 113MeV 

— > rriK 

F T _ >vtt * F r _ >vK' 

Flavor symmetry implies the following relations for the formfactors at t = 0: 

Ff\0) = —F A (0) = -0.0410, 

F^ K) {0) = —F v (0) = -0.0955. (103) 

As far as the quantum numbers are concerned, the resonances if* (892), if* (1410) and if* (1680) 
could contribute in the vector channel and if y (1270) and if \ (1400) in the axial one. The rates 
r(ifi(1400) — ► Kp,Ku) and r(if*(1410) — ► Kp), however, are compatible with zero ||], and so, 
assuming vector meson dominance, their contribution to r — * isKj may be neglected. In the case of 
the pionic mode, the contribution of the p" was found to be small. Therefore we assume that in the 
kaonic mode the contribution of the if*(1680) is small also and approximate it by zero. So we are 
left with the if* (892) and the if i (1270) and a monopole form for both form factors: 

Ff\t) = F ( A K) (0)Re SKl{im (t) 

F^ K \t) = 4 X) (0)Res^ (892) (t). (104) 



For simplicity we use constant widths as in equations ([i5|) and ( |69| ) for both the if* and the ifx- 

The resulting integrated decay rate Tk(Eq) is plotted in Fig. 8 in dependence on the photon 
energy cut Eq. Note that here the structure dependent effects are much more important than in the 
pionic case. This is understood easily, because the higher mass of the kaon suppresses the internal 
bremsstrahlung. So the measurement of structure dependent effects is fairly easy is this decay mode. 
We have also performed the integration using a kaon-photon invariant mass cut. The results are: 

T K (^ = 800 MeV) = 3.58 • 10' 3 T T ^ uK 
T K {^= 1200 MeV) = 0.91 • lO -3 ^-^ (105) 



7 Summary and Discussion 

We have calculated the radiative one pseudoscalar decay modes of the tau lepton. Because of the 
infrared divergence of the internal bremsstrahlung the integrated decay rates depend on a photon 
energy cut Eq or equivalently on a meson-photon invariant mass cut to. The branching ratio for the 
pionic mode is found to be of the order of 10~ 3 (10~ 4 ), using a rather low (high) value for Eq. For 
the kaonic mode, the branching ratio is of the order of 10 -5 . So at least the pionic radiative mode 
should be observable presently. 
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The rate for r — ► vn~{ is dominated by the internal bremsstrahlung, and so, in order to see 
structure dependent effects, sufficient statistics and/or investigation of the invariant mass spectrum 
rather than the photon energy spectrum are needed. Because of the higher mass of the kaon, the 
observation of "structure" will be easier in the decay r — ► uKj. 

Our basic assumption in the parametrization of the form factors Fy and Fa was that they are 
dominated by hadron resonances with suitable quantum numbers. 

In order to parametrize Fy in the pionic case, we then used the CVC-hypothesis to relate Fy(t) 
to i ? 7r7 (t) and the QCD-theorem on lim (tF ni ^ and assumed a saturation of the QCD-theorem 

by the p-like resonances of lowest mass. While the other assumptions seem rather safe, there is no 
obvious a priori justification for the last one. We find, however, that the p(770) alone saturates 
the QCD-theorem at the 86% level and that the inclusion of the p' and the p" changes the picture 
only little, which supports our assumption a posteriori. Then if this saturation is believed in, the 
remaining ambiguities in the parametrization of Fy(t) have only little impact on the decay rate for 
r — > VTTJ. 

In the axial channel the a\ is the only known resonance which has the correct quantum numbers. 
Still FA.(t) dominates the uncertainty of the prediction for r — ► vir^/, because neither Fa(0) from 
radiative pion decay nor the width of the a\ are known very well experimentally. In this paper we 
took the attitude that the negative sign solution for Fa (corresponding to 7 « 1/2) is the physical 
one. Measurements from pions at rest also allow for a solution with 7 ~ —2.5. We showed that the 
prediction for r — > vir^/ changes very significantly if 7 rs —2.5, and so a measurement of this decay 
will allow for a clear independent check of 7 « 1/2. Furthermore we showed that by measuring the 
pion-photon invariant mass spectrum in the a\ region a decision for a small or large a\ width can 
be made. So it will be very interesting to see whether the measurement of r — ► v-K^f will support the 
large r oi (of 400 . . . 600 MeV) as measured in the channel r — > v?>ir or a small a\ width (of about 
300 MeV), as prefered by hadronically produced a\. 

So we believe that once a high statistics sample of r — ► virj is available we will learn much about 
low energy physics. 

In a forthcoming paper we will address the question of soft photons and present the total decay 
rate T(r — ► v T + tt + (7)). 
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